Direct measurement of single soft lipid nano-tubes: nano-scale information extracted 

in non-invasive manner 
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We investigated the dynamics of single soft nano-tubes of phospholipids to extract nano-scale 
information such as size of tube, which are several tens to hundreds of nano-meters thick. The 
scaling law of the nano-tube dynamics in quasi- 2-dimensional space was revealed to be constituent 
with that of a polymer. The dynamic properties of the tubes obtained from direct observation by 
fluorescent microscopy, such as their persistence length, enable us to access the nano-scale char- 
acteristics through a simple elastic model of the membrane. The present methodology should be 
applicable to the nano-sized membrane structure in living cells. 

PACS numbers: 87.16.dj, 87.14.Cc, 87.15.H-, 65.80.-g 
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Within organisms, the phospholipid membrane is a 
primitive compartment wall of small structures that can 
exhibit various topologies and morphologies, such as 
sphere, tube, etc. While there are too many examples 
to enumerate, several characteristic morphologies have 
been reported in endoplasmic reticulum [if, neurons and 
mitochondria [2] . Although the shapes of the membrane 
as observed in electron micrograph images are fixed or 
frozen, they must be dynamic in living cells. For exam- 
ple, the deformation of mitochondria has been observed 
depending on changes in the environment Such dy- 
namic morphologies under non-equilibrium and reactive 
conditions could have biological significance and should 
be examined by physico-chemical methods. This is one of 
the motivations for identifing and measuring the dynam- 
ics of the morphology of nanoscopic membrane structures 
such as tubes. 

Recent biophysical approaches that use model cell sys- 
tems have been shown to be useful for examining the 
above problems. For instance, giant vesicles which ex- 
hibit several tens of micrometers in size have been used as 
a model system for studying their morphological trans- 
formation W, microdomain formation or lateral phase 
separation [5|. The encapsuplation of biopolymers has 
been applied for mimicing biochemical functions 0i9 
and structures 0, ■ Toxicity by an enzyme on a mem- 
brane has been observed as a deformation of the vesi- 



cle [Tl| , and affinity of the DNA on a membrane has been 
identified with transcriptional processes 12, Com- 
plex system or membrane containing biomolecules have 
been realized in supported membrane 14| . 

In the present paper, we report a novel observation 
method for identifying the dynamics of single soft nano- 
tubes, and suggest a new strategy for measuring their 
physical properties or nanoscopic structures in a nonin- 
vasive manner. Dynamical properties, including the dif- 
fusion constant, relaxation time and persistence length, 
are obtained directly through the use of fluorescent mi- 
croscopy. Comparison of a model of membrane elasticity 
based on the Helfrich scheme tlFj to the statics of a fluc- 



tuated tube enables us to determine the thickness of in- 
dividual nano-tubes. This strategy should be applicable 
to non-immersive in vivo measurements within cellular 
organisms. 

Phospholipid tubes were made by the natu- 
ral swelling method [lj|- DOPC (Dioleoyl-L-a- 
phosphatidylcholine, purchased from Wako) and 
Rhodamine-Red DPPE (N-(rhodamine red-X)-l,2- 
dipalmitoyl L-a-phosphatidylethanolamine, Avanti 
Polar Lipids) at a molar ratio of 500 : 1 were dissolved in 
chloroform/methanol (=2:1) solution, and then dried 



(a) 



0.0 



5.0 



10.0 sec 



(A 



5 (.im 




FIG. 1. (a) Typical images of a single phospholipid tube in 
an aqueous solution. The images were taken at intervals of 
5.0 s. The scale bar is 5 /J,m. (b) A schematic illustration of 
a tube and the measured values obtained from image analy- 
sis. The radius of gyration R g is calculated from the spatial 
distribution of fluorescent intensity; the rotational relaxation 
time t> represents the decay of the correlation of a unit end- 
to-end vector p(t); the persistence length l p is defined as the 
decaying slope of a scalar product of two unit tangent vectors 
u(s) and u(s + As), (c) Frame representation of part of a 
tube. An Eel-like tube is assumed to be composed of several 
sections of tori with various values of R and homogeneous r. 
Notations used in the calculations are represented as above. 
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FIG. 2. (a) Log-log plot of measurements of MSD of the 
center of mass. In the short time region to 10 sec, each MSD 
is proportional to At. (b) Log-log plot of MSD of the end 
points. 



under vacuum overnight to form thin lipid films. The 
dry films were hydrated with 200 /iL of distilled water 
at room temperature for 3 hours. A glass chamber for 
microscopy observation was made using glass cover slips 
washed with alkaline alcohol. The swelled solution was 
placed between the slips so that the solution was 2 /im 
thick, and the chamber was sealed with hydrophobic 
liquid blocker to prevent evaporation. The sample 
chamber was set on a fluorescent microscope (Nikon 
TE2000U). Since the sample was confined to within this 
thin space that was on the order of the focal depth, the 
shape and motion of the whole tube could be observed 
as in-focus images. Sequential images were captured by 
an EM-CCD (Hamamatsu Photonics) and recorded at 
31.93 frames/sec. 

Figure QJa) shows typical images of the time devel- 
opment of a single phospholipid tube exhibiting transla- 
tional and intrachain Brownian motion. As shown, we 
observed individual tubes that were longer than several 
micrometers, and measured their positions and shape. 
The center of mass rc.M. of the phospholipid nano-tube 
is calculated as r c .M. = J2i r i ' H r i) / I2i 1 ( r i) , where 
J(rj) is the fluorescent intensity of each pixel r^. 

The mean square displacement (MSD) of rc.M. for each 
tube is shown in Fig. [U(a) as a function of the lag time 
At. From the slopes of MSD in the short time region less 
than 10 sec, the diffusion coefficient Dq is calculated as 
MSD = 4D G At + A(At) 2 , where A is a drift coefficient. 

Figure rjfb) shows the MSD of the end points of each 
tube and the average of scaling indices in the short time 
region less than 10 sec, that is determined to be 0.76 ± 
0.18. The lateral motion of an end point can be predicted 
by a simple scaling law that the index is calculated to be 
2/3 before the configuration undergoes relaxation 




FIG. 3. Relations between (a) the lateral diffusion constant 
Da and the radius of gyration R g and (b) the rotational relax- 
ation time TV and R g as determined experimentally. Fitting 
lines are Dg oc R g a (a = —1.3) and r r oc Rf (/3 = 2.9), re- 
spectively. Superimposed lines show the scaling law from the 
theory of polymer dynamics; a = — 2 and /3 = 4 are indices 
of the Rouse model, and a = — 1 and f3 — 3 are those of the 
Zimm model. 



ment with the model of polymer dynamics. 

As shown in Fig. [ljb), the radius of gyration and 
the unit end-to-end vector are also obtained as Rg = 

EiHrd ■ In - rc.M.lVEiJfa) and P(*) = K=l(*) - 
r s=Q {t))/ \r s=L (t) - r s=Q (t)\, respectively. D G for each 
tube is plotted as a function of R g in Fig. Eta). The 
scaling index a of Dq cx R g a is a = — 1.3 ± 0.2. The 
superimposed lines in Fig. [3](a), a — —2 and a = — 1, 
indicate the indices calculated from the Rouse model and 
Zimm model for polymers, respectively [l7j . The rota- 
tional relaxation time r r of the unit of end-to-end vector 
pit) is also obtained in Fig. UHb) as a function of R g . 
The scaling index f3 for Dq cx R J* in the experiment is 
estimated to be j3 = 2.9 ± 0.2. The indices for the Rouse 
(3 = 4 and Zimm (3 = 3 models are also shown [IT]. The 
indices from the experiment are rather close to those of 
Zimm. 

The experiments show that the quasi-two-dimensional 
chamber enables the observation of the dynamics of 
whole soft nano-tubes. Since, the yielded dyna mics ar e 
similar or comparable to the previous studies 
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and thus the experimental results show quite good agree- 



mm, 

thus the result itself must be reliable. In other words, the 
experimental results show that phospholipid nano-tubes 
behave as a polymer chain of the order of micrometers 
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FIG. 4. Typical examples of the correlations of scalar product 
(it(s) ■ u(s + As)}. The persistence length l p is determined 
from the exponential decay in a short region. The correlation 
at long range and in the vicinity of tends to deviate from 
exponential fitting. 



long. 

Next, we discuss and figure out a simple methodol- 
ogy or strategy for determining the unseeable microscopic 
features of the membrane from observable values by cou- 
pling polymer physics and membrane physics through 
fluctuation- response. 

Geometrically, it is easy to assume a bent tube to be 
comprised of parts of a torus as shown in Fig. [He). The 
bending free energy of the membrane F can be denoted 
by a Helfrich-type expression as 
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F=^J(C 1 +C 2 - CofdA. 



(1) 



For part of a torus, C\ and Ci are curvatures that are 
normal to each other on the membrane surface. Cq is the 
spontaneous curvature and in this case is assigned a value 
of zero. When the small radius and its angle are defined 
as r and 4> for the C\ coordinate, and the large radius 
and its angle are R and for the Ci coordinate, C\ = - 
and Co — p l cos - t are given. The area element dA is 

^ R+r cos tp & 

represented as dA = rd(j> ■ (R + rcos<fi)dQ. Thus, the 
bending energy of a torus for which the cross-sectional 
circle sweeps in small angle AO becomes, 
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where c is the curvature of the circular arc drawn by 
the backbone of the torus, and As = RAO corresponds 
to the small section of the arc-length 
Taylor series of AF torU s for c becomes, 
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The first term on the right-hand side represents the bend- 
ing energy for making a tube from a plane membrane, 
and the second and later terms are the result of bending 
of the backbone. 
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FIG. 5. Distribution of the radii of tubes with k = 20 ksT. 



On the other hand, the bending free energy of a chain 
such as a polymer embedded in two spatial dimensions 
can be represented as 0, Q , 



AF 
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Asc' 



(4) 



where l p is the persistence length of a chain in two spa- 
tial dimensions. If we compare the c 2 terms of eqs. ^ 
and 43} , 
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(•5) 



is obtained as a first approximation. 

The persistence length l p of a nano-tube can be ob- 
tained from the experiments as 



(u(s) ■ u(s + As)) = exp I — 



As 

1 1, 



(6) 



where u(s) is a two-dimensional unit tangent vector at 
which s is the path length along the backbone from one 
end, as shown in Fig. [lib). The value of l p for each tube 
can be estimated from the slope in Fig. [4] Short-range 
correlation lower than 0.5 /im is excluded from the fitting 
analysis [25 , 26 1 . The above equation includes an approx- 
imation which adapts the present geometry and the or- 
der of the persistence length: l p calculated from the two- 
dimensional projection image of the polymer chain within 
a quasi-two-dimensional space asymptotically converges 
to the actual two-dimensional persistence length [27jj . 
The present experimental values, thickness of the confin- 
ing space and the order of magnitude of the persistence 
length, satisfy the above approximation. 

The bending rigidity of a DOPC bilayer membrane has 
been experimentally measured by various methods, such 
as the pipette injection method H, 29} and the electrode- 
formation method (3f3 | , and is known to be approximately 
20 U-qT. Therefore, the thickness of each nano-tube can 
be calculated from eq. (O by substituting the value into 
k, and the result is shown in Fig. [SJ The radii of nano- 
tubes, which are evaluated individually, are distributed 
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around 30 nm. This value is consistent with the mi- 
croscopy observation, considering the limitations due to 
resolution and diffraction limit, where the thickness is 
less than or equal to the fluorescent blurring length. In 
this manner, dynamical analysis of the micrometer-scale 
behaviors has revealed the size of a nano-structure. 

The radii of the tubes can be estimated through an- 
other approach by considering the tube as a polymer 
chain divided into N = L/b segments, where L is the 
contour length and b = 2l p is the Kuhn length. The dif- 
fusion constant of a rod with length b and radius r in two 
spatial dimensions, Ds, is calculated as (l7j 



D s = 



3 fc B rin(6/2r) 

8 TTT] s b 



(7) 



where rj s is the viscosity of bulk water. Ds can also 
be represented as Ds — NDq — LDa/b in the Rouse 
model. Although the compensations for the entire length 
of the tube and the increase in viscosity from the walls 
are considerable 3l|, [32j , a simple expression in terms of 
r can be obtained as 



l p ( 8wrf s 

r \3kbJ 



(8) 



Since only L and Dq are indefinite but experimentally 
observable variables, both r and the bending modulus n 
of the membrane can be estimated from eqs. ([5]) and l[5]). 
However, the calculated values of r obtained from eq. ([5]) 
were scattered due to the exponential term, and the form 
is too sensitive to evaluate r and k here. 

The present discussion enables us to extract informa- 
tion regarding nanometer-scale properties of individual 
tubes from the dynamics, such as the size and bending 
modulus. Under optical microscopy, diffraction limita- 
tions or fluorescent blurring of approximately 200 nm in 
length hides nanoscopic structures. Based solely on a 
spatial analysis, the apparent sizes of the tubes in the 
fluorescent images were close to around 500 nm in diam- 
eter [3^| . On the other hand, the present method can 
access a scale that is one order of magnitude smaller 
than the length of optical resolution. This would be 
advantageous in measurements of living cells compared 
to methods that involve freezing or fixing, e.g., elec- 
tron microscopy, which have thus far been used to mea- 
sure nano-sized structures [33l 34 1. Living cells and 



nonequilibrium conditions also require non-immersive ap- 
proaches. Optical tweezers are useful tools for measure- 
ments in livin g ce lls c omp ared to a tiny needle or mi- 
cropipette 28|, |29|, |35j, |36|, but they can be difficult to 
use with an internal structure confined in a cell. How- 
ever, the present approach is free from such problems 
because of avoiding adhesive manipulations and modifi- 
cations. Thus, this approach enables the observation of 
both living cells and non-equilibrium vesicles. For exam- 
ple, F-BAR protein families which are expected to induce 



budding or tube-generating transition in a living cell also 
make a spherical vesicle into tubules in vitro model ex- 
periment. By use of the present methodology, protein as- 
sociated membrane tubes have been analyzed and yielded 
a mean bending modulus of the membrane reinforced by 
the proteins as a time-development behavior [37| . This 
non-immersive measurement would be a powerful method 
for evaluating nanoscopic sizes and properties in transient 
or nonequilibrium processes. 

Although the model described here is constructed un- 
der the approximation of an elastic membrane, future 
work should consider the implications of the fact that a 
bilayer membrane is a liquid membrane. First, there is 
friction between the membrane and solution. This fric- 
tion is discussed in plane geometry [38[ , but internal fric- 
tion within the tube is more difficult to consider due to 
coupling of the shape of the tube vesicle. Another im- 
portant consideration is thickness uniformity. The bent 
tube in this analysis is assumed to have a uniform thick- 
ness that does not depend on the backbone curvature c. 
However, the bending energy is also a function of the 
tube radius r, and not only c, as in the bending form 
of eq. ([2]) or ([3]). r has a strong effect in higher-order 
terms. Nonlinear elastic behavior of this type is some- 
times observed in experiments, and it would be interest- 
ing to quantify. These factors may affect those on a larger 
scale, i.e., those in eqs. © and ([7]). The nonlinearity of 
single chains might complement the rheology on a macro- 
scopic scale, not only for tubes but also for a worm-like 
micelle solution [39j]. 

Zimm model takes into account hydrodynamic interac- 
tion by considering Oseen Tensor. The relative strength 
of the interaction depends on the thickness compared 
with the persistence length of the tube. Therefore, rela- 
tive hydrodynamic interaction decays between thin seg- 
ments in effectively swelled chains. This would be one of 
reasons of the behavior closing to that of Rouse model. 
Of course, hydrodynamic interaction through the inner 
part of the tube should be effective to the dynamics, but 
it is now an open question. 

In conclusion, this is the first report on the use of direct 
observation to analyze the dynamics of single phospho- 
lipid soft nano-tubes. We have described a method for 
estimating the thickness and rigidity of a phospholipid 
nano-tube from the observed Brownian motion. The 
ability to avoid damaging the sample object should be 
useful not only for measuring the dynamical response of 
a membrane under non-equilibrium conditions but also 
for revealing the internal properties of membrane-closed 
organelles such as mitochondria and endoplasmic reticu- 
lum. 
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DERIVATION OF FREE ENERGY OF 
MEMBRANE 



The non-invasive measurement for nano-meter sized 
objects is extensively studied in nanotechnology and con- 
densed matter physics. Less explored but relevant one 
is membrane nano-tube that is widely seen in synthetic 
lipid vesicles and plasme membranes in living cells. In 
this studt we have developed novel method to infer the 
cross-section of fluctuating single membrane nano-tube. 
This supplemental material describes the brief summary 
for the calculation of the bending energy in eq. (1) shown 
in the manuscript. For general case, spontaneous curva- 
ture Co is regarded as non-zero variable. 

As shown in the manuscript, two principal curvatures 
on the membrane surface are C\ = l/r and C2 = j?™ c s Q t ^ 1 
where r and </> are the polar coordinates representing the 
cross-section of membrane nano-tube and R is the radius 
of the backbone. Then, Eq. (1) of a torus which the 
cross-sectional circle sweeps in small angle AO is written 
as 



unit circle. The residues are given by 

z 2 + 1 

Res (zq) = lim 1 



z-hso [ z(z — Z\)(z — Z2) 
= 1 

f z 2 + 1 
Res (zx) = lim <^ (z - zi) ■ — - 

z^Zl { Z(Z - Zl)(Z - Zi) 

_ zf + l 

z\{zi - z 2 ) 
a 



(S4) 



(S5) 



We use the residue theorem and then rewrite Eq. (|S3[) as 

/ a 2 

■ 2m y Res(z. ( ) = 2tt ( — 



— • 2tti Res (zj) = 2n 
j=o,i 
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Accordingly, we arrive at the following expression for 



eq. UST): 
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a + cos 4> 
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K 2tt 
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(a + cos <j>) , 



(S7) 
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where a = R/r > 1 and b = 1 — Cor. We then rewrite 
Eq. 



AF t , 



kA9/2 



2 71 



(6+1) cos + ab 2 - 



a + cos 4> 



(S2) 

The integrals of first and second terms in Eq. (|S2[) are 
and 2irab 2 respectively. With the complex integration on 
the unit circle C, we obtain the integral of third term as 



{z + z- 1 )/2 dz 
c a + (z + z~ x )/2 iz 
a f z 2 + 1 



i Jc z(z + a - \/a 2 - l)(z + a + V 'a 2 - 1) 



-dz. 



(S3) 



This integrand has two single poles at zq = and z\ = 
— a + \J a 2 — 1 while we do not take into account zi — 
—a — \J a 2 — 1 is not the pole because it is out from the 



If Co = 0, one gets Eq. (2) in the manuscript. We note 
that Cor (Cor — 2) may have no effect in the bending of 
the backbone, so one can apply Eq. (5) written in the 
manuscript even though the spontaneous curvature of the 
membrane Co is non-zero. It just displaces zero point of 
the free energy. 



2. CONFIGURATION AL FLUCTUATION IN 
THE DIRECTION OF DEPTH 

We next discuss additional effects of three-dimensional 
fluctuation of tube and viscosity compensation on the 
present geometry of membrane tube. A possible error 
in inferring the size of membrane tube may come from 
the fluctuation in the direction of depth because it makes 
the contour length L2 that was observed in experiment 
shorter than the actual overall length L. For simplicity, 
a single wave approximation is applied in the present 
geometry as 



z{s) 



Az . ,2-k , 

sin — s). 

2 A 



(0 < s < L 2 ) 



(S8) 
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This displacement is introduced independent of bending 
in the xy-plane. The amplitude 4n and wavelength A of 
displacement are determined by the following consider- 
ations. Since the entire image of a tube is clearly seen 
during observation, the effective thickness is determined 
to be Az = 1 /im. In terms of wavelength A, in nature, 
persistence length l p may correspond to the combination 
of statistical fluctuation in the xy-plane and sine-wave 
displacement z(s). For simplicity one may take a simple 
case with -| = 6. From these considerations, L can be 
represented as 

where Az — 1 /im. In this expression, L is 1.0 to 1.1 
times longer than L 2 obtained from the projected image. 
Another effect is the compensation of viscosity between 



two nearby walls. This can be described for a first ap- 
proximation as 

*(r) = .j^" , (S9) 

1 16 \J ~T d-z> 

, where 77,50, d and z are the viscosity of bulk water, the 
distance between two flat walls and the average distance 
to the central axis of the tube from the wall, respec- 
tively [l|. If we incorporate eq. (|S8|) into rj s of eq. (8) in 
the manuscript, the equation becomes nonlinear equation 
about r. One can get r as well as r) s (r) which becomes 
1.0 to 1.1 times larger than t] s q. 
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